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ABSTRACT 

A  computer-controlled  data  collection  system  was  designed  and 
implemented  to  analyze  finite  amplitude  acoustic  standing  waves  in  a 
rectangular  air-filled  cavity  with  linear  and  wedge-shaped  boundary 
perturbations.  The  response  of  a  nonlinear  standing  wave  in  per- 
turbed and  unperturbed  cavities  was  studied.  The  experimental 
results  showed  that  the  computer-obtained  data  was  in  qualitative 
agreement  with  that  previously  obtained  by  manual  procedures.  Quan- 
titative comparisons  with  theory  could  not  be  made  because  the  prop- 
erties of  the  cavity  did  not  match  those  for  which  theoretical 
predictions  are  available. 
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I.    INTRODUCTION 

A  one-dimensional  mathematical  model  for  finite-amplitude 
standing  waves  in  a  real  cavity  was  studied  by  Coppens  and  Sanders 
[Refs.  1,  2]  in  1968.  In  1974,  this  model  was  extended  to  three 
dimensions.  Losses  and  resonance  frequencies,  inadequately 
described  by  the  classical  theory  of  wall  losses,  are  obtained  for  the 
real  cavity  from  measurements  obtained  in  the  linear  acoustic  regime. 
The  research  done  in  1974  showed  that  when  a  cavity  is  driven  near 
the  resonance  of  a  standing  wave  the  nonlinear  distortion  primarily 
involves  only  the  normal  modes  of  the  cavity  that  are  "familially 
related"  to  the  driven  mode.  Lane  [Ref.  3],  DeVall  [Ref.  41,  Slocum 
[Ref.5],  and  Kuntsal  [Ref.  6]  examined  the  problem  and  showed  that 
experimental  results  are  in  excellent  agreement  with  the  theory 
except  when  a  degeneracy  exists  in  the  cavity. 

The  purpose  of  the  research  reported  in  this  thesis  was  to  design 
and  implement  a  computer-assisted  data  collection  and  analysis  sys- 
tem for  measuring  the  finite-amplitude  distortion  in  a  perturbed 
cavity  with  degeneracies. 
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II.    BACKGROUND  AND  THEORY 

When  the  sound  intensity  in  an  acoustic  process  becomes  so  high 
that  the  principle  of  superposition  no  longer  applies,  we  enter  the 
domain  of  nonlinear  acoustics.  The  limits  of  linear  approximations  in 
the  equation  of  force,  mass  continuity,  and  state  are  exceeded  and  the 
resultant  vibrational  disturbance  at  any  point  in  the  medium  will  no 
longer  be  equal  to  the  sum  of  its  individual  components.  The  essence 
of  the  nonlinear  acoustic  mechanism  is  that  an  intense  wave  perturbs 
the  medium  in  which  it  exists  and  this  perturbation  alters  the  natural 
order  governing  the  wave's  behavior.  The  result  is  that  segments  of 
the  disturbance  in  a  warmer  compression  travel  faster  than  those  in  a 
cooler  rarefaction.  When  this  happens,  the  acoustical  wave  alters  its 
shape,  steepening  as  it  travels  along.  The  alteration  is  infinitesimal  for 
weak  disturbances  and  gradually  increases  in  importance  with 
increases  in  both  amplitude  and  frequency  of  sound.  [Ref.  2] 

In  1968,  Kirchoff  [Ref.  7]  studied  the  distortion  of  intense  acous- 
tic waves  in  an  unbounded  fluid.  This  study  was  extended  by  Lamb 
[Ref.  8],  Fay  [Ref.  9],  and  by  Keller  [Ref.  10]. 

Keck  and  Beyer  [Ref.  11]  approached  the  problem  of  unbounded 
propagation  by  using  perturbation  methods  and  Fourier  series 
representations  of  the  waveform.  The  model  of  Coppens  and  Sanders 
[Ref.  12]  dealt  with  finite-amplitude  standing  waves  in  rigid  walled 
cavities.  They  extended  the  Keek-Beyer  perturbation  approach  to 
include  wall  losses  predicted  by  the  Rayleigh-Kirchoff  theory  [Ref.  7], 
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and  showed  excellent  agreement  with  experiments  conducted  in  a 

rigid-walled    tube    at    low    but    finite    levels    of    excitation.       The 

experimental  work  by  Beech  [Ref.   13]  and  by  Ruff  [Ref.   14]  showed 

that  at  high  excitation  levels  differences  developed  between  theory 

and  experiment.    Experimental  work  by  DeVall  [Ref.  4]  and  by  Kilmer 

[Ref.    15]   showed   that  if  degeneracies   existed   the  model  failed   to 

account  for  the  experimentally  observed  excitation  of  "nonfamilial" 

modes. 

The  effects  of  various  geometrical  boundary  perturbations  on  a 

rigid-walled  cavity  investigated  by  Aydin  [Ref.    16]   showed  that  the 

standing  waves  that  exist  in  an  ideal  cavity  must  be  corrected  when 

the  boundaries  are  irregular.    The  model  of  Coppens  and  Sanders  is 

based    on    a   three-dimensional,    nonlinear   wave    equation   with   a 

dissipative  term  describing  absorptive  losses  encountered  by  standing 

waves  in  rigid  walled  cavities.    As  stated  in  [Ref.  12], 

The  standing  waves  which  can  be  excited  in  any  real  cavity  deviate 
from  predictions  of  the  linear  wave  equation  with  ideal  boundary 
conditions  for  the  following  reasons: 

a     The  presence  of  boundary-layer  losses  at  the  cavity  surfaces 
yields  a  dispersive  contribution  to  the  wave  equation. 

b.     Geometrical  irregularities  alter  the  effective  dimensions  of  the 
cavity. 

Both  of  these  mechanisms  can  be  treated  separately  as  long  as  the 

shift  in  frequency  is  so  small  that  the  actual  resonances  are  close  to 

the  theoretical  values  resulting  from  the  classical  model.    These  are 

treated  by  assuming  the  dimensions  are  exact,  and  the  apparent  phase 

speeds  are  determined  on  that  basis. 
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The  resonance  frequency  for  each  standing  wave  is  defined  as 
[Ref.  6] 


con  =  Cn  [(nx  kx)2  +  (ny  ky)2  +  (nz  kz)2p 


(1) 


nrc 


n\K 


\k 


where  kx  =  r~,    ky  =  -r— ,    kz  =  r—,    and  n,  m,  and  1  are  integers.    Cn  is 

the  apparent  phase  speed  appropriate  for  that  frequency. 

The  force  equation,  in  Eulerian  coordinates,  can  be  written  in  the 
form 


^+fff.vm+ivp  =  ^* 


(2) 


where  the  operator  X  describes  those  physical  processes  leading  to 
absorption  and  dispersion.  If  we  ignore  rotational  effects  so  that 
"u  =  V<t>,  restrict  waves  to  rectangular  cavities,  and  neglect  terms  of 
orders  higher  than  M2,  M(a/k),  and  (cc/k)2,  the  result  is  a 
quadratically  nonlinear  wave  equation  [Ref.  12]. 


p2  -2 
S>  -L 


P    ^ 


PoC^ 


a2 
at2" 


y-  1  (   P 
2     [^C| 


^2 


u 

Co 


(3) 


The  standing  pressure  wave  in  a  rectangular  cavity  has  the  form 
[Ref.  12] 


p  =  coskx  x  coskyy  coskzz  coscot 


(4) 
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If  the  cavity  is  being  driven  near  resonance,  the  nonresonant  terms 
contribute  almost  nothing  with  respect  to  the  resonant  terms.  So  the 
standing  wave  pressure  has  the  form  of 


P  =  X    Pn 
n=l 


(5) 


where  each  Pn  is  of  frequency  nco  and  has  the  form  [Ref.  12]. 


Pn 


: — -pjr  -  MRn  cosnkxx  cosnkyy  cosnkzz  sin(ncot  +  Qn) 


(6) 


These  pns  constitute  those  normal  modes  familially  related  to  pi. 

The  nonlinear,  coupled,  transcendental  equations  which  must  be 
solved  for  the  Rns  are 


Rn 


cos 
sin 


(Qn  -  6n)  =  NM(3Qn  COsQn 


^X    Rj  Rn-j  J  stof  *  (Qj  -  Qn-j) 


COS 

sin 


cos 


I    Rn+jR)lsin[*^n-j-Qj) 


for  all  n  >  1 


(7: 


R=  1        for  n=l 


where  Qn  =  Arctan(-Fn) 


Fn  =  Q  [(nco)2  -  (con)2]/(nco)2  =  2Qn 


C0n 


for 


nco  -  o)n 
con 


«  1     (8) 


9n  = 


C0n 


0>nu  "  WnL 


(9) 


Total  pressure  near  resonance  for  perturbed  boundaries  is  [Ref.  6] 
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p  =  (n,  m,  l/co.0)  - 


(— ? 

C0) 


Z    ami  An  (-Dn  * 


m=o 

1  =  0 


Snml   (n.m,l/co,0  +  Cnml) 


10) 


where    Snmi  -    first-order   perturbation   correction   and   anml    angle 
which  is  defined  as  [Ref.  6]: 


1/fti 


ml 


e  = 


Lx 


A  =  Amp.  of  perturbation 


If  the  cavity  has  perturbed  walls,  the  solution  will  be  in  terms  of  a 
summation  of  the  classical  linear  solution  for  ideal  boundaries  plus 
perturbation  correction  terms  due  to  the  irregular  boundary: 


p  =  p0  +  G  p  +e^p    + 


(11) 


Since  the  magnitude  of  the  boundary  perturbation  is  kept  small,  sec- 
ond and  higher  terms  in  e  can  be  considered  insignificant,  so  that 


p  =  p0  +  e  p 
and  p  must  satisfy  the  following  conditions, 

-lP  =  0 


(12) 


(13) 
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and 


Vp  ■  n  =  0 


at     x  =  0,  Lx  [1  +  ef(y,z)] 

Y    =    0.    Ly 

z  =  0,  Lz 


(14) 


where  n\  the  local  normal  to  the  real  surface,  is  obtained  by  taking  the 
gradient  of  the  equation  for  the  boundary  [Ref.  1], 


n  =  V  {x-Lx  [1  +  ef(y.z)]} 


(15) 


Thus,  to  the  first  order  in  e , 


_^     A       T    3f(y,z)  A       T    3f(y,z) 
n  =x-ELx-^y-€Lx  -£- 


(16) 


and  when  equation  (16)  is  used  in  equation  (15),  the  result  is 


"9£  9f(y,z)  dp  9f(y,z)  dp 

ax_eLx   ay    ay"GLy   az    az. 


x  =  Lx[l+6f(y,z)] 


=  0 


(17) 


Applying  a  Taylor-series  expansion  for  p  evaluated  at  the  real 
boundary  Lx  [1  +  ef(y,z)]  produces 

pi  =  P  I         +  £  I         G  Lxfly.z) 

x  =  Lx[l+EfIy,z)l        x  =  Lx      axx  =  Lx 


+  5-|£  I         [eLxf(y,z)]2  + 

^  CA     x  =  Lx 


(18) 
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Substituting  equation  (12)  into  the  RHS  of  equation  (18),  taking 
the  partial  derivative  with  respect  to  x  on  both  sides  and  keeping  the 
first-order  terms  in  e,  yields  [Ref.  1] 

axx  =  Lx(l  +  efly,z)]      dx  x  =  Lx  dxx  =  Lx       a*    x  =  L*  V 


Taking  the  partial  derivatives  with  respect  to  y  and  z  and  using 
exactly  the  same  procedure  gives  [Ref.  1] 

i£l  -l£l      +*¥"      +^l      «My^)  +  ...    (20) 

^yx  =  Lx[i +  eny.z)i   ^yx  =  Lx      ayx  =  Lx    ^y^xx  =  Lx       u 


^£  I  3P°  I  dp   I  ^2Po  I  T    f ,      v  ro , , 

"3T  =3T  +  €  "37"  +^>7  e  Lxf(y,z)  +  . . .     (21 

^zx  =  Lx[l  +efly.z)l      dzx=Lx  ^zx  =  Lx      ^z^xx  =  Lx  ^ 


Substituting  equations  (19),  (20),  and  (21)  into  equation  (17)  and 
keeping  the  first-order  terms  in  e  results  in 


axx'=Lx-I«L"ily*z)^r+   ay    ay+    az     azjx=Lx 


(22) 


The  RHS  of  equation  (22)  can  be  represented  as  a  Fourier  series  in 
cosines,  so  that  p'  can  be  expressed  as  a  summation  of  normal  modes. 
Hence, 

"■£-1         =   X   I  amicos-r — y  cos  7—  z  cos  (cot  +  0)  (23) 

dxx  =  Lx      m=01=0  ^  ^ 


or 
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TfiH    T     =   II   ami  (O.m,l/w.0)  (24) 

oxx  =  Lx       m=o  1=0 


where 


Ly    L2 

aoo  =  o^J   j   [G]dydz  (2  5; 


2    YY 

amo  =  O^  J    J    [G]  cos  -j—  y  dy  dz 


i_/y   i-2 

J    \   [GJcos^ydydz  (26) 


2      YY  Ltc 

a0l  =  O^J    J    [G]cos^dydz  (27) 


4     V  f   ._.         mic  Ltc  ,     m*0  ,«o* 

a™1  =  L^J  J  'GlcosT7ycosE7zdyd-n^o  (28» 


and  G  is  defined  as 

G.^l  (29) 

In  order  to  find  the  contribution  to  p'  from  each  of  the  terms,  it  is 
stated  [Ref.  16]  that,  for  a  cavity  forced  by  a  dynamic  boundary  condi- 
tion at  a  boundary, 

C^D2Lp'  =  0  (30) 

and  the  dynamic  boundary  condition  from  the  (m,l)th  term  is 
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9^-l_      =  Acos-j— - y  cos-j—z  eu®t+Qj  (3D 


then 


1  °°  T 

I  =  "      ToTS ^    An  l-l)n  Snml  COS  -j^-x  cos^—- y  cos  ^- 


~a\  <* 


.  ei(cot+Q+anmD  (32) 

where 

fl  if  n  =  0 
An  "  [2  if  n=  1,2,3,... 

Applying  this  solution  to  equation  (24)  gives  the  complete  solution 
for  the  first-order  perturbation  [Ref.  6] 


p  =  (n.m.l/co.e)  -77r% —  X   ami  An  (-l)n  Snml 

n=o 
m=o 

1=0 


Co/1^  m=°o 


(n.m.l/w.e+Onml)  (33) 

As  it  is  stated  [Ref.  16],  if  this  is  near  a  resonance,  co-conml.  then 
this  term  will  dominate  the  summation  and  all  other  non-degenerate 
terms  can  be  omitted.    Consequently,  equation  (33)  becomes 

p  =  (n,m,l/co.e)  -777^ — ami  An  (-l)n  Snml  (n,m,l/co,0+anml)  •••    (34) 


-0^ 
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III.    EXPERIMENTAL   CONSIDERATIONS 

A.      COMPUTER-CONTROLLED  APPARATUS 

The  rectangular  cavity  used  in  this  research  was  designed  and 
constructed  in  1977  by  Kuntsal  [Ref.  6]  (Figure  1).  One  wall  of  the 
cavity  could  be  moved  and  secured.  The  interior  of  the  cavity  was 
12.002  inches  long  and  2.502  inches  high,  while  the  width  could  be 
varied  from  approximately  5.50  inches  to  7.00  inches,  in  0.25-inch 
steps.  The  piston,  located  in  a  2.25-inch  diameter  port  in  the  floor  of 
the  cavity,  was  supported  by  an  O-ring  to  produce  a  seal  between  the 
piston  and  the  port.  Three  1/4-inch  diameter  microphone  ports 
were  used.  A  condenser  microphone  was  mounted  in  a  0.24-inch 
aluminum  case  and  an  O-ring  and  silicon  grease  were  used  between 
this  case  and  the  ports  for  sealing. 

Modal  designation  for  the  lowest  three  modes  is  shown  in 
Figure  2.  The  (020)  and  (100)  modes  are  nearly  degenerate.  The 
(020)  mode  is  a  one-dimensional  standing  wave  with  nodes  3  and  9 
inches  along  the  long  wall  of  the  cavity  and  the  (100)  mode  is  a 
standing  wave  with  one  node  in  the  middle  of  the  short  wall.  When 
the  microphone  is  placed  at  position  A,  the  (100)  mode  is  sensed 
without  contamination  from  the  (020)  mode.  In  position  C,  the  (020) 
mode  is  sensed  without  contamination  from  the  (100)  mode.  A  block 
diagram  of  the  apparatus  is  shown  in  Figure  3. 

The  piston  was  driven  by  an  M-B  Electronics  model  EA1500 
exciter,  which  in  turn  was  driven  by  an  APS  Electronic  model   114 
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(a)  Cavity  Orientation 

(b)  Modal  Designation 
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Block  Diagram 
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power  amplifier.  The  driving  signal  was  produced  by  an  HP  3325A 
Synthesizer/Function  generator. 

The  piston  movement  was  continuously  analyzed  by  an  Endevco 
model  2215  Accelerometer  mounted  within  the  piston.  The  output 
voltage  of  the  microphone  and  output  voltage  of  the  accelerometer 
were  measured  by  an  HP  3678  multimeter.  An  HP  3562A  Dynamic 
Signal  analyzer  was  used  to  measure  the  harmonic  distortion.  A  model 
5204  Lock-in  Analyzer  was  used  to  measure  the  resonance  frequen- 
cies at  small  signal  amplitude. 

The  sound  pressure  in  the  cavity  was  sensed  by  a  Bruel  and  Kjaer 
4136  condenser  microphone  with  matching  preamplifier  2807.  The 
output  of  this  preamplifier  went  to  three  pieces  of  equipment:  (1)  an 
HP  3478  multimeter  to  measure  the  voltage  level,  (2)  an  HP  3 5 62 A 
Dynamic  Signal  Analyzer  to  measure  the  harmonics,  and  (3)  a  Cos 
5060A  oscilloscope  to  monitor  the  undesired  distortion  of  the  accel- 
erometer signal. 

The  HP  3562A  Dynamic  Signal  Analyzer  and  HP  3325A  Synthe- 
sizer/Function generator  were  controlled  by  an  HP  87  computer.  Two 
types  of  geometrical  perturbations  were  used: 

a   The   long  wall  was   machined   at   a   small   angle   (as   shown   in 
Figure  4). 

b.  The  long  wall  was  machined  as  shown  in  Figure  5. 

The  effect    of   these    geometrical    perturbations    on    the    resonance 

frequency  of  the  (100)  mode  could  be  estimated  by  calculating  an 

"effective"  width  Le  =  L  -  AL  of  the  perturbed  cavity  where  AL  =  V/A 
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Figure  4 
Linearly  Perturbed  Wall 


.25" 


Figure  5 
Wedged  Perturbation 

and  V  and  A    are    the    volume    and    cross-sectional    area    of   the 
perturbation.  [Ref.  1] 

B.      STRENGTH  PARAMETER  (STRM) 

Strength  of  the  finite  amplitude  interaction  is  characterized  by 
the  strength  parameter  STRM  [Ref.  6]. 

STRM  =  M[3Qi  where  Qi  is  the  quality  factor  of  the  driven  mode 


P  =  2 


U1+Ctf 


M  = 


Pl  = 


V 
Pi 


Cv 


=  1.2  for  air 


pC2 
vi  V2" 

Sm 


where  Sm  is  the  microphone  sensitivity. 
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The  microphone  output  voltage  vi  is  read  on  an  HP  3562A 
Dynamic  Signal  Analyzer  set  at  the  fundamental  harmonic.  Micro- 
phone sensitivity  was  obtained  with  a  Bruel  and  Kjaer  model  4220 
piston  phone. 

SM  =  (1.56  ±  0.03)  x  10-3  volt/(N/M2) 
By  using  p  =  1.293  kg/m3  and  C  =  345  m/sec,  we  found 

STRM  =  7.07  x  10-3  x  V  x  Q.  (35) 

This  is  the  same  sensitivity  as  is  measured  by  Kuntsal  [Ref.  6]. 

C.  FREQUENCY  PARAMETER  (FP) 

The  frequency  parameter  is  defined  by 

FP  =  2Qi  (cd-o>i)/©i  (36) 

where  f  is  a  driving  resonance  frequency  and  fi  is  a  resonance  fre- 

Acon 
quency  of  (010)  mode.     It  is  assumed  that  — —  «   1.     FP  indicates 

where  the  driving  frequency  is  with  respect  to  the  fundamental  reso- 
nance frequency  of  the  cavity.  As  I  FP  I  increases,  the  nonlinear  effects 
will  in  general  become  smaller  for  a  given  strength  parameter. 

D.  HARMONICITY  COEFFICIENT   (E(n)) 

The  harmonicity  coefficient  is  defined  as  [Ref.  6] 

EM  =  ^T^1  (37) 
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E(n)  characterizes  how  well  the  modes  of  a  given  family  are  tuned. 

E.      COMPUTER  PROGRAMS 

Two  computer  programs  were  developed.  The  first  program  finds 
the  resonance  frequencies  and  calculates  the  quality  factors  (Figure  6). 
The  second  program  drives  the  Dynamic  Signal  Analyzer  HP  3562A  by 
a  HP  3325A  Synthesizer/Function  generator  and  finds  and  calculates 
harmonic  content  of  the  microphone  output  at  each  driving  frequency 
(Figure  7). 
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'9"  Calculation  Flow  Diagram 
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Power  of  the  Harmonic  Content  Calculation  Flow  Diagram 
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IV.    DATA  COLLECTION  PROCEDURES 

A.  SYSTEM  ALIGNMENT 

The  primary  factor  in  setting  the  lower  limit  to  the  measurement 
of  harmonic  content  was  the  amount  of  distortion  injected  into  the 
cavity  by  the  piston.  Misalignment  of  the  piston  within  its  port  causes 
distortion  of  the  piston  motion  which  was  detected  by  analyzing  the 
harmonic  content  of  the  accelerometer  output  voltage.  (With  the  pis- 
ton removed  from  its  port,  distortion  levels  were  of  the  same  order  of 
magnitude  as  the  general  electronic  noise  level  of  the  accelerometer 
system  when  at  rest.)  Before  each  run,  the  piston  was  aligned  within 
its  port  so  that  the  accelerometer  output  voltage  contained  less  than 
0.1  percent  second  harmonic  and  0.2  percent  third  harmonic. 

After  testing  various  configurations,  the  one  selected  for  ease  of 
alignment  consisted  of  the  cavity  mounted  on  a  tripod  support  directly 
above  the  exciter,  which  was  supported  by  its  own  base  in  an  upright 
position.  Both  the  cavity  tripod  and  the  exciter  rested  on  a  5/16-inch 
pad  of  stiff  acoustic  rubber.  Once  the  system  was  aligned,  only  an 
occasional  minor  adjustment  was  necessary  to  maintain  low  piston 
distortion. 

B.  AUTOMATIC  DATA  COLLECTION 

Since  the  resonance  frequencies  of  the  cavity  were  observed  to 
drift  with  time,  the  system  was  allowed  to  warm  up  for  at  least  three 
hours  prior  to  data  collection.  The  harmonic  content  of  the 
accelerometer  output  was  analyzed  and  the  piston  was  adjusted  until 


29 


the  second  harmonic  of  the  accelerometer  output  was  at  least  50  dB 

below  that  of  the  fundamental  frequency.    Figure  8  shows  the  drift  in  fr 

is  approximately  linear  in  time. 

Data  was  collected  using  two  measurement  procedures  in  the 

computer  programs:     (1)  a  pre-  and  post-run  infinitesimal  amplitude 

measurement;  and  (2)  the  finite  amplitude  measurement. 

1.      Pre-Run  and  Post-Run  Infinitesimal  Amplitude 
Measurements 

During  the  pre-run  and  post-run  infinitesimal  amplitude 
measurement,  the  piston  was  driven  80  mV  (rms).  Because  the  (100) 
and  (020)  modes  were  nearly  degenerate,  measurements  of  the  prop- 
erties of  the  (100)  mode  were  made  with  the  microphone  in  port  C 
while  measurements  of  the  (020)  mode  were  made  with  the  micro- 
phone in  port  B.  The  resonant  frequency  of  the  fundamental  and  the 
Q  factor  first  few  (0,n,0)  in  the  modes  were  determined  from 

fr  =  ^  (38) 


Q=A:  (39) 


and  the  harmonicity  coefficients  were  then  calculated  from 

fn  ~  nf  i 

e»=V  (40) 
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Figure  8 
Drift  in  Frequency 
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where  fu  and  f^  are  the  frequencies  of  the  upper  and  lower  half-power 
points,  and  fn  is  the  resonance  frequency  of  the  nth  member  of  the 
(010)  family.  Figure  9  shows  the  precision  of  frequency 
measurements. 

2.     The  Finite-Amplitude  Measurement 

All  finite-amplitude  measurements  were  made  at  port  B. 
Throughout  this  measurement,  Vi  had  to  be  kept  constant  by  the 
computer  program  (Appendix  C),  which  required  adjustment  of  the 
driving  voltage  applied  to  the  piston.  The  frequency  was  increased  by 
the  computer  in  0.5-Hz  steps  through  ±5  Hz  about  the  resonance 
frequency  of  a  fundamental  mode.  At  each  driving  frequency,  the 
harmonic  content  of  the  microphone  output,  usually  up  to  the  fourth, 
was  measured  with  the  HP  3562A  Dynamic  Signal  analyzer. 

To  determine  the  exact  value  of  the  frequency  parameter,  the 
value  of  fr  at  the  instant  of  the  measurement  of  Vn  was  determined 
from  the  first  computer  program. 
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Figure  9 
Precision  of  Frequency  Measurements 
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V.    RESULTS 

In  this  section,  the  theoretical  predictions  and  experimental 
results  obtained  by  Kuntsal  [Ref.  6]  will  be  qualitatively  compared  to 
the  experimental  results  obtained  by  using  the  computer-controlled 
equipment. 

The  information  on  the  empirical  losses  and  resonance  frequen- 
cies is  contained  in  the  Q's  and  E's.  It  is  important  to  note  that, 
because  Vi  is  always  kept  constant  as  a  function  of  frequency,  the  sec- 
ond harmonic  distortion  is  maximized  for  co  =  002/2.  The  point  where 
this  occurs  for  each  V2/V1  curve  is  indicated  by  an  arrow  with  the 
label  Fo20-   At  that  point,  the  value  of  the  frequency  parameter  is 

F020  =  2  Q010  E020 
Similarly, 

Fo30  =  2  Q010  Eo30 

is  indicated  on  each  graph.  The  arrow  labelled  F100  indicates  the 
position  where  the  nearly  degenerate  (1,0,0)  mode  is  resonant,  and 
the  value  of  F 100  is 

F100  =  2  Q010  El00 

The  theoretical  predictions  for  unperturbed  cavities  made  by 
Kuntsal  [Ref.  61  were  compared  to  the  experimental  results  as  seen  in 
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Figure  A.  It  is  clearly  seen  that  the  theory  and  experiment  are  qualita- 
tively in  agreement. 

At  Figure  9,  fi  and  (2  are  plotted  as  a  function  of  time.  This  figure 
showed  two  things:  (1)  changes  of  resonance  frequencies  of  the  family 
modes  are  proportional  to  the  frequency,  as  would  be  expected  if  the 
long  dimension  of  the  cavity  changes  with  time,  most  likely  due  to 
temperature  changes;  and  (2)  the  fluctuations  of  resonance  frequen- 
cies from  measurement  to  measurement  indicate  a  precision  of  ±  0.02 
Hz  in  determining  the  resonance  frequencies  of  the  lowest  mode.  By 
using  the  value  of  ±  0.02  Hz,  expected  uncertainties  in  the  E's  are 

fn±0.02  -n(fi  ±0.02) 
AEexpected  ~  nfi 

Calculated  AEexpected  and  measured  AE's  between  pre-  and  post-run 
are  tabulated  in  Table  4.  The  difference  of  the  pre-run  and  post-run 
E's  for  the  family  modes  in  the  unperturbed  cavity  are  within  the 
expected  precision.  For  the  100  mode  (non-family  mode),  this  differ- 
ence is  much  greater  than  expected,  indicating  that  the  cavity  may  be 
changing  its  dimensions  unequally  in  all  directions. 

For  a  linear  perturbation  (Figure  B)  and  for  the  wedge  perturba- 
tion (Figure  C),  the  harmonic  content  in  Reference  6  is  in  qualitative 
agreement  with  experiment,  considering  that  lower  Q's  and  different 
E's  were  observed  compared  to  Kuntsal's  measurement  [Ref.  6]. 
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VI.    CONCLUSIONS 

The  experimental  apparatus,  procedures,  and  computer  programs 
are  capable  of  providing  data  sufficiently  precise  to  test  the  prediction 
of  the  theory  in  the  absence  and  presence  of  geometrical  perturbation. 

Because  the  strength  parameters  (STRM)  and  quality  factors  (Q) 
calculated  in  .the  experiment  did  not  match  the  STRM  and  Q  in  Refer- 
ence 6,  quantitative  agreement  could  not  be  tested.  However,  the 
experimental  results  show  that  a  geometrical  perturbation  alters  the 
finite  amplitude  behavior  of  the  cavity,  and  that  the  nature  of  these 
changes  is  in  qualitative  agreement  with  the  predictions  of  the  theory. 

The  quantitative  disagreements  between  the  theory  and  experi- 
ment for  the  perturbed  cavities  may  be  because  the  resonance  fre- 
quency of  the  non-family  mode  is  not  tracking  the  family  modes. 

In  order  to  avoid  disturbing  the  position  of  the  piston  in  the  port 
of  the  cavity,  the  microphone  has  to  be  shifted  very  carefully  from  one 
port  to  another  between  the  "pauses"  as  the  program  goes  on. 
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APPENDIX  A 
CURVES 
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Figure  A 
Unperturbed  Cavity 
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Figure  B 
Linear  Perturbed  Cavity 
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Figure  C 
Wedge  Perturbed  Cavity 
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APPENDIX  B 
TABLES 

TABLE  1 

UNPERTURBED  CAVITY 

(SEE  FIGURE  A) 

Pre -Run 


Mode 

fu 

fL 

fr 

9 

E 

010 
020 
030 
100 

565.801 
1131.840 
1698.271 
1131.872 

562.875 
1127.791 
1693.524 
1127.590 

564.351 
1129.812 
1695.901 
1129.731 

195.40 
278.69 
356.51 
264.20 

0.0 
98.16x10-5 
16. 82xl0-4 
91.17x10-5 

Vi  =  -11  dB 


Finite-Amplitude  Run 
Mode:    010 


STRM  =  .389 


f  (Hz) 

V2  (dB) 

V3  (dB) 

V4  (dB) 

V2/Vi 

V3/V1 

V4/V1 

FP 

558.5 

-50.0 

-70.5 

-70.5 

0.0113 

0.0014 

0.0012 

-4.19 

559.0 

-47.9 

-72.0 

-69.0 

0.0143 

0.0013 

0.0011 

-3.85 

559.5 

-46.9 

-73.8 

-67.8 

0.0160 

- 

0.0011 

-3.50 

560.0 

-46.0 

-72.6 

-67.5 

0.0179 

- 

0.0010 

-3.15 

560.5 

-44.8 

-71.2 

-69.7 

0.0204 

0.0011 

0.0010 

-2.81 

561.0 

-43.7 

-69.7 

-73.5 

0.0233 

0.0012 

- 

-2.47 

561.5 

-42.4 

-67.8 

-76.5 

0.0269 

0.0014 

- 

-2.12 

562.0 

-41.0 

-65.7 

-75.9 

0.0316 

0.0018 

- 

-1.78 

562.5 

-39.4 

-63.2 

-82.5 

0.0380 

0.0025 

- 

-1.43 

563.0 

-37.8 

-60.3 

-85.0 

0.0457 

0.0034 

- 

-1.09 

563.5 

-35.9 

-56.6 

-74.9 

0.0569 

0.0052 

- 

-0.74 

564.0 

-34.0 

-52.2 

-67.5 

0.0708 

0.0070 

0.0015 

-0.40 

564.5 

-32.2 

-47.6 

-60.6 

0.0876 

0.0148 

0.0033 

-0.05 

565.0 

-31.1 

-45.0 

-54.5 

0.0994 

0.0200 

0.0067 

0.29 

565.5 

-31.6 

-43.9 

-52.9 

0.0933 

0.0226 

0.0070 

0.64 

566.0 

-33.3 

-46.7 

-56.0 

0.0772 

0.0164 

0.0056 

0.98 

566.5 

-35.3 

-51.0 

-63.2 

0.0610 

0.0100 

0.0025 

1.33 

567.0 

-37.6 

-55.5 

-70.8 

0.0484 

0.0060 

0.0010 

1.67 

567.5 

-39.8 

-60.0 

-73.6 

0.0365 

0.0036 

- 

2.02 

568.0 

-42.1 

-64.3 

-74.7 

0.0279 

0.0022 

- 

2.36 

568.5 

-44.1 

-68.5 

-74.2 

0.0221 

0.0013 

- 

2.71 

569.0 

-46.2 

-69.5 

- 

0.0174 

0.0011 

- 

3.05 

40 


O 

6" 

12" 
Post-Run 

Mode 

fu 

fL 

fr 

Q 

E 

010 
020 
030 

100 

566.381 
1132.740 
1699.770 
1132.800 

563.270 
1128.676 
1694.974 
1128.487 

564.82 
1130.71 
1697.37 
1130.64 

194.03 
277.82 
353.84 
262.03 

0.0 
94.00xlO-5 
16.90xl0-4 
74.16x10-5 
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TABLE  2 

LINEAR  PERTURBATION  CONFIGURATION 

(SEE  FIGURE  B) 

Pre -Run 


Mode 

fu 

fL 

fr 

9 

E 

010 
020 
030 
100 

569.465 
1133.883 
1703.014 
1131.650 

566.545 
1129.961 
1698.008 
1131.341 

565.653 
1132.306 
1700.594 
1134.651 

193.72 
288.71 
339.72 
263.321 

0.0 
17.67x10-4 
21.42x10-4 
29.57x10-4 

Vi  =  -12  dB 


Finite-Amplitude  Run 
Mode:    010 


STRM  =  .355 


f  (Hz) 

V2  (dB) 

V3  (dB) 

V4  (dB) 

V2/V1 

V3/V1 

V4/V1 

FP 

561.0 

-43.0 

-68.7 

0.0251 

0.0013 

-3.19 

561.5 

-42.7 

-68.0 

- 

0.0260 

0.0014 

- 

-2.84 

562.0 

-42.3 

-66.9 

- 

0.0271 

0.0016 

- 

-2.50 

562.5 

-41.6 

-64.5 

- 

0.0296 

0.0021 

- 

-2.16 

563.0 

-40.3 

-62.7 

- 

0.0342 

0.0026 

- 

-1.82 

563.5 

-38.7 

-60.1 

- 

0.0411 

0.0035 

- 

-1.47 

564.0 

-37.3 

-56.2 

- 

0.0485 

0.0055 

- 

-1.13 

564.5 

-34.2 

-53.5 

- 

0.0691 

0.0075 

- 

-0.79 

565.0 

-33.8 

-51.0 

-65.0 

0.0721 

0.0100 

0.0020 

-0.44 

565.5 

-31.8 

-46.4 

-60.1 

0.0911 

0.0170 

0.0035 

-0.10 

566.0 

-30.9 

-43.0 

-54.6 

0.1012 

0.0251 

0.0066 

0.24 

566.5 

-29.6 

-40.3 

-50.0 

0.1161 

0.0344 

0.0112 

0.58 

567.0 

-30.9 

-45.8 

-54.1 

0.1000 

0.0180 

0.0070 

0.92 

567.5 

-32.8 

-49.8 

-59.2 

0.0810 

0.0115 

0.0039 

1.26 

568.0 

-34.1 

-51.5 

-65.0 

0.0700 

0.0094 

0.0020 

1.61 

568.5 

-36.7 

-56.2 

-68.1 

0.0514 

0.0055 

0.0014 

1.95 

569.0 

-38.0 

-58.7 

- 

0.0441 

0.0041 

- 

2.29 
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Post-Run 


Mode 

fu 

fL 

fr 

9 

E 

010 
020 
030 

100 

570.410 
1135.067 
1705.321 
1136.740 

567.456 
1131.116 
1700.340 
1132.416 

566.861 
1133.984 
1701.110 
1135.756 

191.90 
287.01 
341.50 
262.66 

0.0 
23.10x10-4 
39.98x10-4 
17.94x10-4 
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TABLE  3 

WEDGE  PERTURBATION  CONFIGURATION 

(SEE  FIGURE  C) 

Pre -Run 


Mode 

fu 

fL 

fr 

9 

E 

010 
020 
030 

100 

566.912 
1127.191 
1695.310 
1126.761 

563.864 
1123.281 
1690.471 
1122.911 

565.381 
1125.230 
1692.891 
1124.829 

185.80 
288.08 
349.27 
291.41 

0.0 
-48.92x10-4 
-19.17x10-4 
-52.47x10-4 

Vi  =  -11  dB 


Finite-Amplitude  Run 
Mode:    010 


STRM  =  .370 


f  (Hz) 

V2  (dB) 

V3  (dB) 

V4  (dB) 

V2/V! 

V3/V1 

V4/V1 

FP 

559.4 

-39.5 

-66.9 

0.0376 

0.0016 

-3.99 

559.9 

-37.9 

-65.3 

- 

0.0450 

0.0019 

- 

-3.66 

560.4 

-36.4 

-63.9 

- 

0.0540 

0.0023 

- 

-3.33 

560.9 

-34.5 

-61.8 

- 

0.0672 

0.0029 

- 

-3.00 

561.4 

-32.4 

-59.1 

-78.0 

0.0856 

0.0039 

- 

-2.66 

561.9 

-30.3 

-55.7 

-72.7 

0.1090 

0.0058 

0.0008 

-2.33 

562.4 

-28.7 

-55.6 

-68.2 

0.1310 

0.0059 

0.0014 

-2.00 

562.9 

-28.9 

-50.8 

-65.3 

0.1271 

0.0103 

0.0017 

-1.67 

563.4 

-30.9 

-50.0 

-62.5 

0.1010 

0.0112 

0.0027 

-1.34 

563.9 

-32.3 

-49.7 

-60.3 

0.0810 

0.0116 

0.0034 

-1.01 

564.4 

-34.7 

-50.9 

-61.4 

0.0650 

0.0101 

0.0030 

-0.67 

564.9 

-36.3 

-54.3 

-67.0 

0.0550 

0.0068 

0.0016 

-0.34 

565.4 

-37.5 

-59.3 

-75.2 

0.0480 

0.0038 

- 

-0.01 

565.9 

-38.4 

-65.0 

- 

0.0430 

0.0020 

- 

0.32 

566.4 

-39.3 

-69.7 

- 

0.0391 

0.0012 

- 

0.65 

566.9 

-39.7 

-72.2 

- 

0.0370 

0.0009 

- 

0.98 

567.4 

-40.6 

-74.5 

- 

0.0280 

- 

- 

1.32 

567.9 

-40.7 

-75.4 

- 

0.0230 

- 

- 

1.65 

568.4 

-41.6 

-73.5 

- 

0.0210 

- 

- 

1.98 

568.9 

-44.0 

-71.5 

- 

0.0193 

- 

- 

2.31 

569.4 

-47.3 

-70.9 

- 

0.0151 

- 

- 

2.64 

569.9 

-49.1 

-71.1 

- 

0.0120 

- 

- 

2.97 
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025"' 

e- 

J 

6" 

O 

12" 

Post-Run 

Mode 

fu 

fL 

fr 

9 

E 

010 
020 
030 
100 

566.950 
1127.311 
1695.544 
1126.981 

563.963 
1123.411 
1690.731 
1123.126 

565.464 
1125.371 
1693.144 
1125.051 

189.24 
288.18 
351.57 
291.39 

0.0 
-49.12x10-4 
-27.33x10-4 
-50.21x10-4 
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TABLE  4 
E'S  DIFFERENCES  BETWEEN  RUNS 


Unperturbed 
Configuration 

Linear  Perturbed 
Configuration 

Wedge  Perturbed 
Configuration 

AE2 

.416xl0"4 

5.43x10-4 

0.2x10-4 

AE3 

.08xl0-4 

8.56x10-4 

8.16x10-4 

AEioo 

1.7xl0-4 

11.63x10-4 

2.26x10-4 

AEexpected 

.53x10-4 
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APPENDIX  C 
COMPUTER  PROGRAMS 

This  program  drives  the  HP  3325A  Synthesizer/Function  genera- 
tor between  562  Hz  and  1710  Hz,  finds  the  resonance  frequency,  and 
calculates  the  quality  factor  (Q). 

10  !  SUBROUTINE  Q 

14  SETTIME  0,87200 

1  5  FOR  INC=1  TO  30 

20  K=-1 

30  FOR  FREQ=562  TO  1 71 0  STEP  2 

50  IF  FREQ>550  AND  FREQ<570  THEN  GOTO  210 

70  IF  FREQ>570  THEN  GOTO  90 

90  FOR  FREQ=1 1 1 6  TO  1 1 40  STEP  2 

110  IF  FREQ>1 1 15  AND  FREQ<1 140  THEN  GOTO  210 

1 30  IF  FRQ>1 1 40  THEN  GOTO  1 791 

140  K=K+1 

150  FOR  FREQ=1 680  TO  1 71 0  STEP  2 

1 70  IF  FREQ>1 670  AND  FREQ<1 71 0  THEN  GOTO  21 0 

190  IF  FREQ>1 71 0  THEN  GOTO  1 800 

210  F1=FREQ 

230  F2=F1+2 

250  F3=F1+4 

270  REMOTE  717,723,702 

290  OUTPUT  71 7  ;"FR";F1  ;"HZ" 

310  WAIT  5000 

330  ENTER  723  ;  AMP1 

350  OUTPUT  71 7  ;"FR";F2;"HZ 

370  WAIT  5000 

390  ENTER  723  ;  AMP2 

41  0  OUTPUT  717  ;MFR";F3;"HZ" 

430  WAIT  5000 

450  ENTER  723  ;  AMP3 

470  DISP  F1  ,F2,F3,AMP1  ,AMP2,AMP3,K 

51 0  IF  AMP2<AMP1  OR  AMP2<AMP3  THEN  GOTO  590 

570  IF  AMP2>AMP1  OR  AMP2>AMP3  THEN  GOSUB  BINOM 

590  NEXTFREQ 

610  !  SUBROUTINE  BINOM 
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630      BINOM: 

650      P=F2 

670      S=2 

690      S=S/2 

710      OUTPUT717;"FR";P;"HZ" 

720      WAIT  5000 

730      ENTER  723 ; AP 

750      IF  P>1 71 0  THEN  GOTO  1 800 

770      PU=P+S 

790      OUTPUT  71 7  ;"FR";PU;"HZ" 

800      WAIT  5000 

810      ENTER  723  ;  AU 

830      PL=P-S 

850      OUTPUT  71 7  ;"FR";PL;"HZ" 

860      WAIT  5000 

870      ENTER  723 ; AL 

890      DISP  P,PU,PL,AP,AU,AL 

910      IFAU>APTHENP=PU 

930      IFAL>APTHEN  P=PL 

950      IF  S>.008  THEN  GOTO  690 

970      GOSUB  HPWR 

990      RETURN 

1010    !  SUBROUTINE  HALFPOWER 

1030   HPWR: 

1050    DISP  P,AP 

1060    DISP  TIME 

1065    PRINT  TIME 

1070    PRINT  "MAX" 

1090    PRINT  P.AP 

1100    IF  INC>1  AND  INC<30  THEN  GOTO  1 771  ELSE  GOTO  1110 

1110    PRINT  "REAL  HALFPOWER" 

1130   APH1=AP/2A.5 

1150    PRINT  APH1 

1 1  70    FOR  F=1  TO  2 

1190    PH=P+((-1)AF) 

1210   OUTPUT717;"FR";PH;"HZ" 

1230   WAIT  5000 

1250    ENTER  723; AH 

1270    DISP  PH.AH 

1290    IFAH<APH1  THEN  GOTO  1350  ELSE  GOTO  1310 

1310    PH=PH+((-1)AF) 

1330   GOTO  1210 

1350   S=.5 

1370    PHI=PH-(S*((-1)AF)) 

1390    FOR  1=1  TO  10 

1410   OUTPUT717;"FR";PHI;"HZ" 

1430   WAIT  5000 
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1450  ENTER  723  ;AHI 

1470  S=S/2 

1490  IF  AHI>APH1  THEN  GOTO  1570 

1510  PHI=PHI-(S*((-1)AF)) 

1530  NEXT  I 

1550  GOTO  1610 

1570  PHI=PHI+(S*((-1)AF)) 

1590  NEXT  I 

1610  PRINT  "HALFPOWER" 

1630  DISP  PHI.AHI 

1650  PRINT  PHI.AHI 

1670  IF  F=1  THENPHIR=PHI 

1960  NEXTF 

1710  Q=P/(PHI-PHIR) 

1730  PRINT  "Q  VALUE" 

1750  PRINT  Q 

1770  DISPQ 

1771  K=K+1 

1772  IF  K=0  THEN  GOTO  90 

1  773  !  AFTER  FIRST  PAUSE  SHIFT  THE  MICROPHONE  TO  THE  PORT  C  TO 

MEASURE  THE  FREQUENCY  OF  THE  DEGENERATE  MODE. 

1  775  !  AFTER  SECOND  PAUSE  SHIFT  THE  MICROPHONE  TO  THE  PORT  B 

1791  IF  K=1  THEN  PAUSE 

1792  IFK=2THEN  PAUSE 

1793  IF  K=1  THEN  GOTO  90 

1794  IF  K=2  THEN  GOTO  140 
1797  RETURN 

1800  NEXT  INC 

1810  END 


49 


This  program  drives  the  HP  3562A  Dynamic  SignalAnalyzer  by  HP 
3325A  Synthesizer/Function  generator  and  finds  and  computes  the 
harmonic  content  of  the  microphone  output  at  each  driving  frequency 
between  563  Hz  and  570  Hz. 

10  REMOTE  720,717 

15  SETTIME  0,87033 

20  OUTPUT  720  ;"PRST" 

30  OUTPUT  720  ;"B" 

40  OUTPUT  720  ;"FRS3KHZ" 

50  OUTPUT  720  ;"PSPC" 

60  P=800 

70  OUTPUT  71 7  ;"AMP";P;"MV" 

80  OUTPUT  720  ;"LNRS" 

100  OUTPUT  720  ;"STBL" 

1 1 0  OUTPUT  720  ;"PSP2" 

1  20  OUTPUT  720  ;"KEY31  " 

1  25  !  DRIVING  THE  DYNAMIC  SIGNAL  ANALYZER  AND  READING  THE 

HARMONIC  CONTENT  OF  THE  DRIVING  FREQUENCY. 

130  FREQ-563 

1  40  OUTPUT  71 7  ;"FR";FREQ;"HZ" 

1  50  OUTPUT  720  ;"  SF563HZ  " 

160  WAIT  3000 

1  70  OUTPUT  720  ;"STRT" 

1  80  WAIT  8000 

1  90  OUTPUT  720  ;"MKPK" 

200  OUTPUT  720  ;"RDMK" 

210  ENTER  720  ;X1,YB 

220  DIL=YB 

230  DISPX1.DIL 

234  PRINT  TIME 

235  PRINT  X1.DIL 

240  FOR  FREQ=563  TO  570  STEP  .5 

250  OUTPUT  720  ;HFRS3KHZ" 

251  OUTPUT  720  ;"SF";FREQ;"HZ" 
255  OUTPUT  71 7  ;"FR";FREQ;"HZ" 
260  FOR  1=1  TO  4 

270  OUTPUT  71 7  ;"AMP";P;"MV" 

280  OUTPUT  720  ;"STRT" 

290  WAIT  8000 

300  OUTPUT  720  ;"MKPK" 
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310  OUTPUT  720  ;"RDMK" 

320  ENTER  720  ;  X,YB 

324  DISPX.YB 

325  PRINT  X.YB 

326  PRINT  TIME 

332  IF  1=1  THEN  GOTO  333  ELSE  GOTO  350 

333  IF  FREQ>563  THEN  GOTO  335  ELSE  GOTO  350 
335  GOSUBBINOM 

350  OUTPUT  720  ;"SF";X+50;"HZ" 

351  OUTPUT  720  ;"  FRS.7KHZ" 

355  NEXT  I 

356  NEXTFREQ 

360  !  SUBROUTINE  BINOM 

365  !  THIS  SUBROUTINE  KEEPS  THE  FUNDAMENTAL  CONSTANT. 

370  BINOM: 

380  RAT=1 0A(DIL/1 0)/1 0A(YB/1 0) 

385  DISP  RAT.DIL.YB 

390  IF  RAT>.999  AND  RAT<1 .001  THEN  RETURN 

400  P=P*RATA.2 

41  0  OUTPUT  71 7  ;"AMP";P;"MV" 

420  OUTPUT  720  ;"STRT" 

430  WAIT  8000 

440  OUTPUT  720  ;"MKPK" 

450  OUTPUT  720  ;"RDMK" 

460  ENTER  720  ;  X.YB 

470  DISPX.YB 

475  PRINTTIME 

480  PRINT  X,YB 

490  IF  FREQ>570  GOTO  660 

510  GOTO  380 

660  END 


51 


LIST  OF  REFERENCES 

1.  Coppens,  A.  B.,  private  communication. 

2.  Sanders,  J.  V.,  private  communication. 

3.  Lane.  C.  Finite  Amplitude  Waves  in  a  Rigid  Walled  Cavity,  Thesis. 
Naval  Postgraduate  School.  Monterey,  California,  1972. 

4.  DeVall,  R.  R.,  Finite  Amplitude  Waves  in  Imperfect  Cavities,  Thesis. 
Naval  Postgraduate  School.  Monterey,  California,  1973. 

5.  Slocum  IV,  W.  S.,  Finite  Amplitude  Standing  Waves  in  Real  Cavities 
Containing  Degenerate  Modes,  Thesis,  Naval  Postgraduate  School, 
Monterey,  California.  1975. 

6.  Kuntsal,  E.,  Finite  Amplitude  Standing  Waves  in  Rectangular  Cavi- 
ties with  Perturbed  Boundaries,  Thesis,  Naval  Postgraduate  School, 
Monterey,  California.  December  1978.    . 

7.  Kirchoff.  G..  Ann.  Phys.  Leipzig,  Vol.  134.  pp.  177-193.  1868. 

8.  Lamb,  H.,  Dynamical  Theory  of  Sound,  2nd  ed..  Ch.  6,  Edward 
Arnold  and  Co.,  London.  England.  1925. 

9.  Fay,  R.  D.,  "Successful  Method  of  Attack  on  Progressive  Finite 
Waves."  J.  Acoust.  Soc.  Am.,  Vol.  28.  pp.  910-914,  September 
1956. 

10.  Keller,  J.  B.,  "Finite  Amplitude  Sound  Produced  by  a  Piston  in  a 
Closed  Tube,"  J.  Acoust.  Soc.  Am.,  Vol.  26,  pp.  253-254,  1954. 

11.  Keck.  W.,  and  R.  T.  Beyer,  "Frequency  Spectrum  of  Finite  Ampli- 
tude Ultrasonic  Waves  in  Liquids,"  Phy.  Fluids,  Vol.  3,  pp.  346- 
352,   1960. 

12.  Coppens,  A.  B..  and  J.  V.  Sanders,  "Finite  Amplitude  Standing 
Waves  within  Real  Cavities."  J.  Acoust.  Soc.  Am.,  58  (6),  December 
1975. 

13.  Beech,  W.  L..  Finite  Amplitude  Standing  Waves  in  Rigid  Walled 
Tubes,  Thesis,  Naval  Postgraduate  School,  Monterey,  California, 
1967. 


52 


14.  Ruff.  P.  G.,  Finite  Amplitude  Standing  Waves  in  Rigid  Walled  Tubes, 
Thesis.  Naval  Postgraduate  School,  Monterey,  California.  1967. 

15.  Kilmer,  M.  J.  Ill,  Finite  Amplitude  Effects  in  Rectangular  Cavities 
with  Perturbed  Boundaries,  Thesis,  Naval  Postgraduate  School, 
Monterey,  California,  1975. 

16.  Aydin.  M..  Theoretical  Study  of  Finite  Amplitude  Standing  Waves  in 
Rectangular  Cavities  with  Perturbed  Boundaries,  Thesis,  Naval  Post- 
graduate School,  Monterey,  California,  1978. 


53 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  VA    22304-6145 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  CA    93943-5002 

3.  Deniz  Kuwetleri  Komutanligi  5 
Bakanliklar,  Ankara/Turkey 

4.  Deniz  Harp  Okulu  Komutanligi  Kutuphanesi  1 
Tuzla,  Istanbul/Turkey 

5.  S.  L.  Garrett,  Code  61GX  1 
Underwater  Acoustics 

Naval  Postgraduate  School 
Monterey,  CA    93943-5000 

6.  Professor  J.  V.  Sanders,  Code  61Sd  1 
Physics  Department 

Naval  Postgraduate  School 
Monterey,  CA    93943-5000 

7.  Professor  A.  B.  Coppens,  Code  6 lCz  1 
Physics  Department 

Naval  Postgraduate  School 
Monterey,  CA    93943-5000 

8.  LT  R.  D.  Middleton,  Jr.  1 
141  John  Street 

Hauppauge,  NY    11787 

9.  Ltjg  Hakan  Basaran  3 
Cemil  Topuzlu  Cad.  Cinar  Apt. 

No.  15/11 

Fenerbahce,  Istanbul/Turkey 

10.  Chief  of  Naval  Operations  1 
(Antisubmarine  Warfare) 

Code  OP-02 

Washington  D.C.   20350-2000 

54 


Thesis 

B2423465  Basaran 

c#1        c°mputer-controlled 
measurements  of  non- 
linear standing  waves. 


H  KOV  90 


3  6  0  9  P 


Thesis 

B2423465  Basaran 

c#l        Co^.putp.r-controlled 
measurements  of  non- 
linear standing  waves, 


NB|«2 


£«A  «-£ 


